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Long Gamma-Ray Bursts (GRBs) are the most dramatic examples of massive stellar deaths, usually as- 
sociated with supernovae [1]. They release ultra-relativistic jets producing non-thermal emission through 
synchrotron radiation as they interact with the surrounding medium [2]. Here we report observations of 
the peculiar GRB 101225A (the “Christinas burst”)* Its 7-ray emission was exceptionally long and followed 
by a bright X-ray transient with a hot thermal component and an unusual optical counterpart During the 
first 10 days, the optical emission evolved as an expanding, cooling blackbody after which an additional 
component, consistent with a faint supernova, emerged. We determine its distance to 1.6 Gpc by fitting the 
spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep 
optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a helium 
star-neutron star merger that underwent a common envelope phase expelling its hydrogen envelope. The 
resulting explosion created a GRB-Iike jet which gets thermalized by interacting with the dense, previously 
ejected material and thus creating the observed black-body, until finally the emission from the supernova 
dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy [3]. 

On Dec. 25, 2010, 18:37:45 UT (T 0 ), the Swift Burst Alert Telescope (BAT, 15-350 keV) detected GRB 101 225 A, 
one of the longest GRBs ever observed by Swift [4] with T 90 > 2000 s (the time in which 90% of the 7 -ray energy 
is released, [5]). A bright X-ray afterglow was detected for two days and a counterpart in ultraviolet, optical and 
infrared bands could be observed from 0.38 h to two months after the event (see the Supplementary Information, 
SI). No counterpart was detected at radio frequencies [ 6 , 7], 

The most surprising feature of GRB 101225 A is the spectral energy distribution (SED) of its afterglow. 
The X-ray SED is best modeled with a combination of an absorbed power-law and a black-body (BB). The 
UV/optical/NIR (UVOIR) SED (see Fig. 1) can be modeled with a cooling and expanding BB model until 10 
days post burst, after which we observe an additional spectral component accompanied by a flattening of the light 
curve (Fig. 2). This behaviour differs from a normal GRB where the SED follows a power-law due to synchrotron 
emission created in shocks when the jet hits the interstellar medium (e.g., [ 2 ]). 

An optical spectrum taken two nights after the burst does not show any spectral lines (see SI). Were fit the SED 
and light curve with the template of SN 1998bw, a SN Type Ic associated with GRB 980425 [44], and obtain a 
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redshift of z = 0.33. At this distance, the SN has an absolute peak magnitude of only My — —16.7 mag, which 
makes it the faintest SN associated with a long GRB [9, 10]. In contrast, the 7-ray isotropic-equivalent energy 
release at z = 0.33 is > 1.4 x 10 51 erg, typical for other long GRBs but more luminous than most other low- 
redshift GRBs associated with SNe [11]. We detect a possible host galaxy in g’ and r' bands with OSIRIS/GTC at 
6 months after the burst with an absolute magnitude of only M 5 , a b S = — 13.75 mag, ~ 2 mag fainter than any other 
GRB host [12]. Although its blue colour matches that of a star-forming galaxy, our observations do not allow us 
to resolve it as an extended source. 

At z = 0.33, the X-ray BB has a radius of ~ 2 x 10 11 cm (~1 R Q ) and a temperature of ~ 1 keV (10 7 K) 
at 0.07 d with little temporal evolution. Such a thermal component, attributed to the shock breakout from the star, 
had also been observed for XRF 060218 [13], XRF 100316D/SN 2010dh [14] and GRB 090618 [15], all nearby 
GRBs associated to Type Ic SNe [16, 17, 14, 18], with similar temperatures but larger radii. The UVOIR BB starts 
with a radius of 2 x 10 14 cm (~ 13 AU) and a temperature of 8.5 x 10 4 K at similar times and evolves considerably 
over the next 10 days reaching a radius of 7x 10 14 cm and temperature of 5,000 K. The evolution of the two BB 
components suggests that they must stem from different processes and regions. 

An appealing model is a Helium - neutron star (NS) merger with a common envelope (CE) phase, a model 
that had been proposed earlier as a possible progenitor for GRBs [19, 20, 21]. In this scenario, a binary system 
consisting of two massive stars survives the collapse of the more massive component to a neutron star (NS). When 
the second star leaves the main sequence and expands, it engulfs the NS, leading to a CE phase and the ejection of 
the hydrogen envelope and part of the helium core as the remnant spirals into the center of the second star. When 
the NS reaches the center, angular momentum forms a disk around the remnant of the merger, allowing for the 
formation of a GRB-like jet. The remnant might be a magnetar whose prolonged activity can explain the very long 
duration of the GRB. 

The interaction of this ultra-relativistic, well-collimated jet with the previously ejected CE material can explain 
both the X-ray and UVOIR emission components. Estimating that the inspiral takes 5 orbits or 1 .5 yr and material 
is ejected at escape velocity, the outer ejecta are at a few 10 14 cm at the time of the merger, consistent with the 
radius of the UVOIR BB. We assume that the ejecta form a broad torus with a narrow, low -density funnel along the 


3 



rotation axis of the system that permits the passage of the 7 -radiation generated in the jet. Most of the jet hits the 
inner boundary of the CE-ejecta and only a small fraction of it propagates through the funnel. The X-ray emission 
is produced by shocks created by the interaction of the jet with the inner boundary of the CE shell. As the jet 
passes through the funnel, it decelerates due to the increased baryon load and shear with the funnel walls so that 
no regular afterglow signature is produced. When the now mildly relativistic, mass-loaded jet breaks out of the 
CE-ejecta, it produces the UVOIR emission in the first 10 d. As the supernova shock expands beyond the CE shell 
we observe a small bump in the lightcurve at ~30 d. The He-NS merger scenario naturally assumes a relatively 
small Ni-production, leading to a weak SN. 

A somewhat similar scenario might also explain XRF 060218 albeit with a different progenitor system pro- 
ducing a brighter SN and a fainter GRB. The remnant might be a magnetar whose prolonged activity can explain 
the very long duration of the GRB. A thermal component in the optical can also be fitted for the X-ray out- 
burst associated with SN 2008D, a Type lb SN in NGC 2770 [22], which, however, had no thermal component 
in X-rays. And last, it is possible that a GRBs shows a thermal component in X-rays (e.g., GRB 090618) but 
shows a classical, bright power-law shaped afterglow [18]. GRB 101225A might be a member of a newly defined 
class of “blackbody-dominated”, SN-associated long-duration GRBs arising in very dense environments created 
by the progenitor systems themselves, which thermalizes the high-energy output from the collapsing star. The 
non-relativistic, uncollimated emission from this scenario makes it difficult to detect them at higher redshifts. This 
makes GRB 101225A a fortunate coincidence that allows us to derive conclusions about the progenitor system and 
its environment by a new variety of massive stellar death, which had so far been only proposed to exist theoretically. 
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frequency (Hz) 

Figure 1: Temporal evolution of the UVOIR spectral energy distribution. The UV, optical and IR counterpart were detected 
by UVOT, the UV telescope onboard Swift and several ground based facilities from 0.38h to nearly 2 month after the GRB. 
This plot shows the evolution of the SED from the onset of the optical observations at 0.07 days to 40 days for all epochs with 
sufficient data to model the SED shape. Filled circles are detections, triangles mark upper limits. The orange line on top of 
the BB model at 2.0 days shows our flux-calibrated spectrum taken with the OSIRIS/GTC. The SED evolution requires two 
different components, a simple expanding and cooling BB up to 10 days and an additional SN component for the last 3 epochs. 
The solid line shows the combined evolution taking both the BB and the SN beyond 5 days into account, whereas the dashed 
line shows the evolution of the BB component alone. The UVOIR BB evolves from an initial temperature of 43,000 K (0.07 d) 
to 5,000 K (18 d) and increases in radius from 2x 10 14 cm to 7 x 10 14 cm at the same timescale. The SN component is best fit 
with a template of the broad-line Type Ic SN 1998bw which was associated to GRB 980425. For the fit of the SN component 
we used the SED at 40 days. Reanalyzing UVOIR data of XRF 060218 [13] and SN 2008D [22] we find a similar thermal 
component over the first 3-4 days after which the SN starts to dominate. The UVOIR BB of XRT 060218 shows a similar 
radius, temperature and evolution as GRB 101225A and is equally inconsistent with radius and temperature from the X-ray BB. 
SN 2008D shows a steeper evolution in the radius, consistent with the continuation from the shock breakout, a thermal X-ray 
component was not observed (see SI). 
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time since burst (days) 

Figure 2: Light curves of GRB 101225A in X-ray and UV/optical/IR bands, (a) X-ray lightcurve, (b) light curves in UV 
to IR bands. Filled points are detections, triangles upper limits. The solid lines denote the combined light curve from the BB 
and the SN component, excluding the contribution from the host galaxy. The evolution of the BB component is shown as a 
dashed line, the SN as a dotted line. Observations started almost simultaneously in X-rays and optical/UV wavelengths. The 
X-rays reached a peak flux of 4.34 x 10~ 9 erg cm 2 s' 1 . After an initial shallow decay of slope £- 1108 ± 0011 U p to 21 ks, 
the X-rays show a strong decay with a slope of £ -8 ' 95±0 ' 20 , inconsistent with synchrotron emission. The UVOIR light curves 
have a shallow maximum at the beginning, with different peak times for the individual bands due to the maximum of the BB 
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emission passing through the spectrum. The second component emerging at around 10 days post-burst is the contribution of 


an underlying SN, modeled with the GRB-SN 1998bw as a template, stretched in time by a factor of 1.25 and decreased in 
luminosity by a factor of 12 (in restframe). The absolute luminosity of the SN is Mv = —16.7 mag, the faintest SN associated 
with a GRB. At ~ 180 days we detect the very faint host at magnitudes of g = 27.36 ± 0.27 and r' = 26.90 ± 0.14 or 
M a hs — —13.7 (0.001 L : ). (c) Color image of the field of GRB 101225A observed at 40 days (indicated by a grey bar in panel 
b) with the afterglow marked by a box. 
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Supplementary information 


1 BAT data analysis and fitting 

GRB 101225A was detected by BAT on-board the Swift satellite [23] on Dec. 25, 2010 at To = 18:37:45 UT as 
an image trigger [24]. It was already in progress both when the source entered the BAT field of view and when it 
left the field of view due to Swift orbit-constrained slews [25, 26]. Therefore, we can only give lower limits on the 
total burst fluence and the T90 duration. The total fluence of the intervals covered in the observations adds up to 
(5.6 ± 0.7) x 10“ 6 erg cm -2 (implying a total energy release in 7-rays of E 7)iso > 1.4 x 10 51 erg at z = 0.33.), 
which is a lower limit to the total gamma-ray emission. No emission was detected in a previous observation of the 
field at T - T 0 = -4950 s and we therefore put a lower limit on the duration of T90 > 2000 s. This is one of 
the highest durations ever observed for a GRB, comparable to the longest burst observed by Swift, GRB 0904 17B 
[27]. The BAT light curve is shown in Fig. 1. The BAT-observed peak flux of (3.25+0.47) x 10~ 9 erg cm -2 
s~ x in the 15 - 150 keV range occurred in the interval T - To = +1372 to T - T 0 = +1672 s. No other 7- 
ray instrument detected GRB 10 1225 A, although the MAXI instrument on board the ISS (2 — 10 keV) reported a 
marginal detection at T — Tq = +1002 s coincident with the BAT position [28]. 

The time-averaged spectra from To to T — To = +963 s and from T — To = +1372 to T — To = +1672 s 
are best fit by simple power-law models with photon indices of T = 1.9 ± 0.4 and 1.9 + 0.2, respectively. The 
lotal fluences in each time period in the 15 - 150 keV band taking these fits are (1.7 ± 0.4) x 10 -6 and (9.0 + 0.2) 
xlO -7 erg cm~ 2 . All quoted errors are at the 90% confidence level. The BAT spectra are almost equally well 
parametrised by models using a cutoff power-law or a blackbody fit due to the low signal-to-noise ratio of the 
event. E pea k using a cutoff model is poorly constrained to 38 ± 20 keV. The blackbody fit gives a temperature of 
kT = 10.1 ± 1.1 keV. Errors are at the 68% confidence level. 

We also examined the BAT data to search for persistent emission after the trigger. For this we used the daily 
sky image mosaics produced as part of the BAT hard X-ray transient monitor which cover a single energy band 
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Figure 1: BAT 7 -ray lightcurve in the 15-150 keV band. The sum of the flux was obtained using a fixed power- 
law index of F = 1.87 from a fit to the most intense part of the burst. The burst started before the beginning of the 
BAT data at ~ T — T 0 = — 100 s and probably continued while the source was not in the BAT field of view from 
T - T 0 = +1091 to T - T 0 = +1372 s. The latest upper limit before the burst was 2.65 x 10 -9 erg cm -2 s -1 at 
T — T 0 = -4950 s. Error bars are at 90% confidence. 


of 15 — 50 keV. We found a 5.3<r excess (0.0048 ± 0.0009 count cm -2 s -1 ) on Dec. 25, 2010 (MJD 55555), the 
day of the trigger, and a positive excess in the count rate (> la or 0.001 1 count cm -2 s -1 ) over the next ten days 
(until MJD 55565). We determine the probability that such a sequence of excess rates would occur by chance. To 
do this, we examine the light curves of 106 “blank sky” points tracked in the BAT transient monitor. These are 
points chosen randomly across the sky at least 10 arcmin from any known X-ray source. Any positive flux from 
these points is expected to be due to chance fluctuations. In these 106 light curves (> 200, 000 data points), we 
find only one sequence of six consecutive days showing a positive excess and none with more than six days. The 
chance probability of ten days of excess flux is less than 1/200,000, so the observed prolonged emission is likely 
real. 

2 XRT data analysis and fitting 

The Swift XRT data were processed with version 3.7 of the XRT data reduction software (released as part of 
HEASoft 6.10 on 2010-09-28) and the corresponding calibration files used for subsequent spectral analysis. 
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The object was detected by XRT from 1.4 ks to 10 5 s after the trigger. Data were collected in Windowed Timing 
(WT) mode for the first 7.3 ks after the trigger followed by Photon Counting (PC) mode for the rest of the obser- 
vations. The peak flux in X-rays is 4.3 x 10 -9 erg cm -2 s'" 1 , the total observed fluence 8.2 x 10 -6 erg cm -2 and 
the unabsorbed fluence 1.1 x 10 -5 erg cm -2 . At z — 0.33 this corresponds to a total energy release in X-rays of 
3.6 x 10 51 erg. Spectra were extracted for individual snapshots of data (one snapshot corresponds roughly to one 
orbit constrained by the observability of the object during the orbit) and were further timesliced into 1 00 s bins for 
the initial snapshot (1.4 — 1.8 ks after the trigger). 



frequency (Hz) 


Figure 2: Fit to the X-ray spectrum. Panel a): X-ray spectrum from XRT during the first snapshot, the dashed 
line indicates the contribution of the power-law component, the dotted line shows the BB component. Panel b): 
Ratio between the observed data and the fitted model. 

We tried a variety of fits to the X-ray data, using XSPEC version 12.6.0, with the result that an absorbed 
power-law plus BB component provided a good fit to the data (see Fig. 2). For the fit we used the Tiibingen-Boulder 
absorption model with the Wilms abundances [29] and Vemer absorption cross-section [30]. As shown in Fig. 3, 
there is little spectral evolution within the first snapshot of data. The best fit taking the entire first 367 s of data is 
a power-law with a photon index F 1.83lo.io* a blackbody of temperature 0.96 ± 0.13 keV (1.11 x 10 7 K) and 


13 





rest frame time since BAT trigger (s) 


Figure 3: Results from the fits to the first snapshot of the X-ray data. The panels show from top to bottom: a) 
count rate during the first snapshot, b) photon index T, c) BB Temperature in keV, d) the contribution of the BB 
to the total emission in percent, e) radius of the emitting BB at z = 0.33 and f) total absorbing column density in 
X-rays (the Galactic column density in the line-of-sight is 7.9 x 10 20 cm -2 ). 


a total absorbing column of (2.2 ± 0.3) x 10 21 cm -2 , for a x 2 of 420.7 for 379 degrees of freedom. The Galactic 
column density in this direction is 7.9 x 10 20 cm -2 . 

The inclusion of the BB is significant at the > 99.9999% level, the contribution of the BB to the total emission 
is around 20% (see Fig. 3). The second snapshot of data (also in WT mode) is again better fit with a BB in addition 
to the power-law, with F = 2. 181^09’ black-body kT = 0.991 q;J 7 keV and NH = (2.7 ± 0.2) x 10 21 cm -2 , 
with x 2 /d.o.f = 378/421. This BB is significant at 99.987%. For the X-ray data after the second snapshot, no 
BB component is required and a simple absorbed power-law provides an acceptable fit, likely due to the lower 
signal-to-noise ratio at later times. 
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Table 1: Result of different model fits to the X-ray data of the first snapshot. PL is a pure power- law model, 
PL+BB a combination of a simple PL and a BB with one temperature while a PL+diskBB includes emission from 
BBs with different temperatures, PL+compt includes a comptonized component in addition to the power-law. T is 
the photon index of the power-law, kT the BB temperature, the column density is the total density, including the 
Galactic absorption. The last column shows the x 2 of the different fits and the F-test value compared to the simple 
absorbed PL model. 


Model 

r 

kT (keV) 

N h (10 22 cm" 2 ) 

X 2 /d.o.f. (F-test) 

PL 

PL+BB 

PL+diskBB 

PL+compt 

1.72±0.03 

i oT+o.ia 
1.50_g 10 

l-79±g;|| 

1.79±0.04 

0.96±0.13 

1.64±0.35 

<22 

0.24±0.02 

0.22±0.03 

0.18+0.04 

O.Ol_0-02 

468/381 

421/379(1. 95 xl0“ 9 ) 
417/379 (3.19xl0 -1 °) 
424/377 (1.55xl0- 7 ) 


The detection of N# in excess of the Galactic column density can be used to constrain the redshift since 
column density and redshift are coupled. We fitted the spectrum of the first orbit fixing the Galactic absorption 
to 7.93 x 10 20 cm -2 . In addition to the absorbed powerlaw plus BB component as described above we allowed 
for a third component leaving both the column density and the redshift free to vary (in contrast to the fit described 
above where the redshift was fixed to 0.3). The resulting contours are plotted in Fig. 4. The 99% upper limit for 
the redshift is 0.5, the 90% upper limit on the redshift is 0.35, consistent with our findings from the UVOIR SED 
fitting. 

We also checked for a possible periodicity in the X-ray data. To this end, light curves were extracted with 18 ms 
bins which is the best time resolution available for WT mode. Using the Kronos powspec tool, no significant 
periodic signal was identified with a frequency between 0.005 and 28 Hz (0.04 and 200 s) in either the first or 
second snapshot of data. The 3 a limit on the variation in flux of any periodic signal during the first WT snapshot 
is 6.3%. 

3 UY, optical and IR data analysis 

3.1 SwiftflTVOT 

Swtft/UVOT began observing GRB 101 225 A 1373 s after the BAT trigger, simultaneous with the XRT observa- 
tions. The automatic target sequence did not start until the end of the BAT image trigger at ~ 23 minutes. The 
optical counterpart was found to be blue, with strong detections in the UV filters (uvw 1, uvm2, uvw2), weak 
detections in the b and u filters, and no detection in the v filter. The data were processed using the standard Swift 
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Figure 4: Redshift vs. absorbing column density confidence contours. The fit is based on the data from the first 
snapshot. Plotted are the 1,2 and 3 a contours together with the 90% and 99% confidence interval. 

software tool uvotmaghist within HE As of t 6 . 9 and the latest calibration files (20101231). 

We extracted counts using a circular aperture with a radius of 5 arcsec for count rates above 0.5 counts s -1 , 
and a 3 arcsec aperture where the count rate drops below 0.5 counts s _1 , as well as a source-free background 
region. The tool uvotmaghist applies coincidence-loss corrections and aperture corrections. The count rates 
were converted to flux density using standard photometric calibration [31, 32]. 

3.2 McD 2.1m 

The CQUEAN instrument (Camera for QUasars in the EArly uNiverse; Park et al. 2011, in preparation) on the 
2.1m Otto-Struve telescope at McDonald Observatory, Texas, USA, observed the optical counterpart starting at 
01:16:23 UT, on Dec. 26, 2010 (6.64 h after the burst). Three exposures of 300 s each were taken in r', i', z' , 
and Y bands under photometric conditions. The data were reduced with standard procedures of dark and flat-field 
corrections. The afterglow is detected in the r\ i' , and .z'-band images, the Y band gives an upper limit only. 
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Table 2: Log of observations in the UVOER range. Values are not corrected for Galactic extinction, 


Mid Hq 
( days) 

Exposure 

(s) 

Filter 

Telescope 

“•'MB 

Flux 

(MJy) 

Premag ing 

3x500 

s' 

3.5mCFHT 

>26.9 (27.2 ±0.5) 

<0.06 <0.048 ±0.22) 

Premseing 

3X500 

r' 

3.5mCFHT 

>25.5 

:0.22 

0.01848 

168 

w2 

UVOT 

> 21.36 

< 10.38 

0.07041 

1431 

w2 

LTVOT 

2156± 0.20 

8.65 ± 1.80 

0.17373 

6719 

w2 

UVOT 

21,63±0.1l 

8.06± 0.83 

0.30739 

6679 

w2 

UVOT 

21.76± 0.12 

7.19± 0.85 

0.45280 

5805 

w2 

UVOT 

21.96± 0.15 

5.96 ± 0.91 

0.81302 

12039 

w2 

UVOT 

2257 ± 0.17 

3.42 ± 0.58 

1.00869 

11753 

w2 

UVOT 

22.37 ± 0.16 

4.08 ± 0.67 

1.41736 

23440 

w2 

UVOT 

22.61 ± 0.20 

3.27 ± 0.66 

1.75211 

23368 

w2 

UVOT 

23 .45 ±0.30 

1.51 ±0.49 

2.44862. 

74516 

w2 

UVOT 

: 23.73 

< 1.17 

4.07964 

13B747 

w2 

UVOT 

>24.20 

<0.76 

7.52954 

377712 

w2 

UVOT 

: 25.39 

<0.25 

0.01818 

319 

m2 

UVOT 

>20.81 

< 17.14 

0.07515 

1431 

m2 

UVOT 

21.97± 0.31 

5.89 ± 1.94 

0.61452 

897 

m2 

UVOT 

21.90± 0.21 

6.34 ± 1.33 

0.95369 

12104 

m2 

UVOT 

22.00 ± 0.15 

5,74 ± 0.83 

1.18487 

18396 

m2 

UVOT 

22.47 ± 0.23 

3.74 ± 0.87 

1.48360 

23468 

m2 

UVOT 

22.46 ± 0.19 

3.76 ± 073 

1.85629 

29528 

m2 

UVOT 

22.97 ± 0.22 

235 ± 0.53 

2.51507 

40973 

m2 

UVOT 

23.34 ± 0.30 

1.68 ± 0.53 

4.08285 

138701 

m2 

UVOT 

>24.25 

:0.73 

0.01846 

318 

wl 

UVOT 

>21.15 

< 12.64 

0.07752 

1431 

wl 

UVOT 

: 22.10 

<5.26 

0.65205 

5571 

wl 

UVOT 

21.81 ± 0.17 

688 ± 1.15 

0.96984 

16520 

wl 

UVOT 

21.72± 0.16 

7.46 ± 1.19 

1.37145 

18904 

wl 

UVOT 

22.23 ± 0.26 

4.65 ± 1.25 

1.71425 

29052 

v! 

UVOT 

22.46 ± 0.25 

3.76 ± 0.97 

2.44380 

74329 

«I 

UVOT 

>22.97 

< 2.36 

4.07673 

138760 

wl 

UVOT 

>23.22 

< 1.86 

0.01789 

169 

u 

UVOT 

>20.46 

< 23.68 

0.07228 

2579 

u 

UVOT 

21.59± 028 

8.42 ± 247 

1.28374 

103571 

u 

UVOT 

22.33 ± 026 

4.24± 1.14 

2.44533 

74215 

u 

UVOT 

>21.82 

<6.80 

4.07781 

138647 

u 

UVOT 

>22.06 

<5.46 

0.01817 

169 

b 

UVOT 

>19.94 

( 38.37 

0.06802 

1430 

b 

UVOT 

7 20.86 

< 16.51 

0.16319 

6726 

b 

UVOT 

2l.53± 0.30 

8.83 ± 2.86 

0.30647 

8349 

b 

UVOT 

>21.83 

<6.76 

1.14624 

127458 

b 

UVOT 

>22.19 

<4.85 

2.44615 

74224 

b 

UVOT 

>21.00 

< 14.47 

4.07834 

128632 

b 

UVOT 

7 21.33 

< 10.64 

0.01789 

318 

1> 

UVOT 

>19.35 

<66.30 

0.07278 

1431 

1' 

UVOT 

7 70 ,36 

C 26.12 

0.18317 

6538 

l- 

UVOT 

: -21,04 

< 13.96 

0.31290 

5819 

V 

UVOT 

>20.69 

< 19.20 

1.11739 

21X180 

V 

1.23mCAHA 

22.47 ±0.19 

3.73 ±0.65 

1.18728 

121098 

i* 

UVOT 

>21.08 

■; 13.39 

2.45087 

74279 

V 

UVOT 

>20.50 

< 2293 

4.08128 

138537 

V 

UVOT 

2 20,92 

: 15J6 

39.11207 

6x180 

g r 

OSIRlS/104mGTC 

> 26.3 

<0.11 

39.49403 

5X180 

g ( 

GMOS/KinOemuli 

26.80 ±0.35 

0.07 ±0.03 

180 

42:: 200 

s' 

OSIRIS/ 10. 4mOTC 

27.21 ±027 

0.047 ±0.0 10 

1.04545 

19x180 

R 

1.23mCAHA 

22.61 ±0.16 

3.28 ±0.48 

0.29887 

3 -.300 

s 

CQUEAN/2. ImMcD 

22.43 ±0.14 

3.87 ±0.50 

108833 

1X30 

/ 

OSIRIS/ lCUmCTC 

23.39 ±012 

1.60±0.18 

21.15017 

101160 

r' 

OSIRIS/ 10.4mGTC 

24.21 ±0.14 

0.75 ±0.10 

28.49818 

5X180 

r 1 

GMOS/8 mGe mini 

24.81±0.13 

043 ±0.05 

39.10159 

4X120 

/ 

OSIRIS/ 1 0.4mGTC 

24.77 ±0.1 3 

0.45 ±0.05 

39.47981 

5 : 180 

r> 

GMOS/8 mGe mini 

25.24±0.15 

0.29 ±0.04 

44.08258 

4X1C0 

/ 

OSIRIS/ 10.4mGTC 

> 24.7 

<0.48 

~180 

32 X 200 

S 

OSIRIS/ 10.4mGTC 

26.90±0.14 

0.063 ±0.008 

LI 7359 

17x180 

i 

1.23mCAHA 

2218±075 

4.88±1.57 

61.96267 

20x120 

l 

SCORPlO/6mBTA 

25.17±0.35 

0.31 ±0.10 

0.29516 

3x300 

l' 

CQUEAN/2. ImMcD 

2272 ±0.1 8 

2.96±0.49 

10.09449 

9x900 

/' 

RAT/2.0mLT 

24.01±0.I3 

0.90±0.11 

39.12164 

5x60 


OSIRIS/ 1 0.4mGTC 

24.36 ±0.17 

0.65 ±011 

3946326 

5X180 

i' 

GMOS/8 mGe mini 

24.61 ±0.09 

052±004 

0J0384 

3x300 

z' 

CQUEAN/2. ImMcD 

2265 ±0.34 

3.16± 1J00 

39.09432 

6X60 

J 

OSIRIS/ 10.4mGTC 

24.73±0.42 

0.47 ±0.1 8 

39.44619 

7x180 


GMOS/8 mGe mini 

24.77 ±025 

0.45 ±0.10 

0J0745 

3x300 

Y 

CQUEAN/2. ImMcD 

> 225 

: 3.63 

37.45092 

32X60 

J 

NlRI/SmGemini 

>23.4 

< 158 

28.46873 

44x60 

K s 

NIRI/8mGeitrini 

24.48 ±035 

059 ±0.19 


3.3 CAHA 1.23m 

The 1.23m telescope is located at the German-Spanish observatory of Calar Alto (CAHA) in Almeria, Spain and 
is equipped with an optical imaging camera. The optical counterpart was detected in the VRI bands 1.04 — 1.11 
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days after the GRB trigger. The 1 .23m was also used to calibrate the object field in BVRI bands by observing the 
Landolt fields RU149D and SA98 on Dec. 26 and 27, 2010, under photometric conditions. 

3.4 LT 2.0m 

The Liverpool telescope is a 2.0m fully robotic telescope located at the observatory of Roque de los Muchachos 
on La Palma. Observations were carried out with the imaging camera RATCAM. The optical counterpart was 
detected in one epoch in i'-band at 10.09 d after the GRB. 

3.5 OSIRIS/GTC 

We acquired imaging data using OSIRIS at the Gran Telescopio de Canarias (GTC), a 10.4m telescope located at 
the observatory of Roque de los Muchachos on La Palma, Canary Islands, Spain. 

The observations started in r' band ~2 days after the burst, exposing for 30 s. A second r' observation was 
carried out ~ 21 days after the gamma-ray event, where 5 exposures of 180 s were obtained. We furthermore 
obtained a late- time SED at 39 days in </, r' , i' and z' bands and a last image at ~44 days in the r' band. 

The data of our two late epochs (at ~39 and ~44 days) were obtained at a considerable airmass (1.73 - 2.14) 
since the object was setting quickly after evening twilight. The data at ~21 and ~44 days were acquired with the 
Moon at ~54 and ~36 degrees, respectively, with an illumination of 83% and 21%, respectively. The SED at 39 
days was constructed in dark time. The observing conditions were good in all four GTC epochs. 

3.6 Gemini-North: NIRI and GMOS-N 

Late-time imaging of the optical counterpart of GRB 101225 A was conducted with the Gemini-North observatory 
on Mauna Kea/Big Island, Hawaii, on several occasions. On the night of Jan. 23, 2011 we observed the field with 
the Near InfraRed Imager (NIRI) in the K' filter for 44 x 60 s exposures (2 x 30 s co-adds) before switching to 
the Gemini Multi-Object Spectrograph (GMOS-N) for 5 x 180 s exposures each in the r' filter. On the night of 
Feb. 1, 201 1 we re-observed the field with NIRI in the J band for 32 x 60 s exposures (co-added), and finally on 
the night of Feb. 3, 201 1 we imaged the field in all four GMOS broad-band filters {g'r'i'z'). Since the source was 
setting, all exposures were taken at moderate to high airmass (1.5 — 2.5), although under relatively good seeing 
conditions. 
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3.7 BTA 6m 


A final late image was obtained using SCORPIO on the 6.0m BTA telescope, located at the Special Astrophysical 


Observatory, in Russia. The observation consisted of 20 x 120 s exposures using an I filter obtained on Feb. 25, 


201 1 (60 days post burst) under good weather conditions and a seeing of 1.3 — 2.0 arcsec. 


3.8 Photometry of ground-based data 


The photometry of V, R and I data was done in a consistent way using a set of 15 comparison stars in the field 
of GRB 101225A calibrated with the Landolt fields taken on Dec. 26 and 27, 201 1 by the 1.23m CAHA. For /, 
i' and z' photometric calibrations was done with observations form the 2.1m Otto-Struve telescope at McDonald 
Observatory on Dec. 26, 2010, using the standard star Feige 34. Finally, g' photometry was derived from the 
other reference magnitudes using numerical transformations [33], The magnitudes of the comparison stars in the 
different filters used for the optical observations are listed in Tab. 3. 

We performed aperture photometry using PHOT within IRAF taking an aperture radius equal to the Full Width 
at Half Maximum (FWHM) of the stellar point sources. In cases where the contamination by neighboring sources 
was not negligible, we did PSF photometry within IRAF. In Tab. 2 we list the final photometry for all UV, optical 
and IR data. 


Table 3: Magnitudes of calibration stars used for the optical photometry. A finding chart indicating the position of 
each reference star, is given in Fig. 5. V, R and I are given in Vega system, the other bands are in AB. 


ID 

R.A. (J2000) 

Dec. (J2000) 

V 

R 

/ 

g' 

~P 


z r 

1 

00:00:48.13 

+44:38:18.9 

15.282 ± 0.009 

14.916 ± 0.049 

14.553 ± 0.052 

— 

— 

— 

— 

2 

00:00:52.02 

+44:37:46.6 

18.065 ± 0.097 

17.268 ± 0.055 

16.288 ±0.110 

18.62±0.23 

17.679±0.01 1 

17.121±0.010 

16.870±0.010 

3 

00:00:41.58 

+44:37:18.0 

16.794 ± 0.050 

16.421 ± 0.051 

15.972 ±0.062 

17.14±0.12 

— 

— 

— 

4 

00:00:52.13 

+44:37:07.9 

18.789 ± 0.079 

18.335 ±0.133 

18.051 ±0.162 

19.26±0.19 

18.464±0,012 

18.325±0.011 

18.263 ±0.013 

5 

00:00:45.46 

+44:37:06.6 

14.995 ± 0.008 

14.633 ± 0.049 

14.235 ±0.051 

— 



. — 

— 

6 

00:00:46.18 

+44:36:52.4 

15.174 ± 0.015 

14.838 ± 0.049 

14,456 ± 0.051 

— 

— 

— 

— 

7 

00:00:55.37 

+44:36:36.8 

19.165 ±0.180 

18.268 ± 0.090 

17,356 ± 0.188 

— 

18.569±0.012 

17. 905 ±0.0 11 

17.597±0.011 

S 

00:01:03.76 

+44:36:26.8 

16.594 ±0.012 

16.161 ± 0.049 

15.691 ±0.090 

— 

— 

. — 

— 

9 

00:00:48.48 

+44:36:19.3 

18.796 ± 0.059 

17.924 ± 0.069 

16.694 ±0.092 

19.08±0.14 

18.602±0.012 

17,819±0.011 

17.516±0.011 

10 

00:00:47.98 

+44:35:57.8 

18.682 ±0.132 

' 17.901 ±0.107 

17.059 ±0.182 

19.41±0.31 

18.174±0.011 

17.645±0.010 

17.427±0.012 

11 

00:00:50.58 

+44:35:43.5 

18.505 ±0.218 

17.877 ± 0.060 

17.270 ±0.163 

— 

18.088±0.011 

17.715±0.010 

17.537±0.012 

12 

00:00:51.59 

+44:35:19.1 

19.253 ± 0.130 

18.893 ±0.131 

18.460 ± 0.273 

19.54±0.31 

19.059±0.014 

18.910±0.013 

18.852±0.020 

13 

00:00:43.29 

+44:35:13.1 

18.050 ± 0.066 

17.749 ± 0.085 

17.318 ±0.108 

18.34±0.16 

17.853±0.011 

17.729±0.010 

17.671±0.012 

14 

00:01:03.09 

+44:35:04.8 

16.828 ±0.015 

16.321 ± 0.052 

15.826 ± 0.087 

— 

— 

— 

— 

15 

00:00:54.83 

+44:35:01.9 

17.516 ± 0.032 

17.024 ± 0.066 

16.321 ± 0.082 

18.04±0.08 

17.148±0.010 

I6.960±0.010 

16.853±0.011 
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Figure 5: Field with the secondary standards used for the photometric calibration. The position of the optical 
counterpart of GRB 101225A is indicated with an arrow. The photometric magnitudes of each of the reference 
stars are given in Table 3. The field of view is 6 arcmin x 4 arcmin. 


3.9 Preimaging 

Preimaging of the field was available from the archive of the 3.5m Canada-France-Hawaii Telescope (CFHT), 
obtained with the MegaPrime/MegaCam for the Pan- Andromeda Archaeological Survey (PAndAS, [34]). We 
combined 3 x 500 s exposures obtained under very good conditions in g' and i' bands and derive 3cr limiting 
magnitudes for these exposures of i' > 25.5 and g' > 26.9. At the position of GRB 101225A, we tentatively 
detect an object with g' — 27.2 ± 0.5 (2er, see Fig. 6) which is consistent with our late host detection (see Sect. 
3.10). 

3.10 Host detection 


At three epochs, June 9, June 28 and 29, we obtained deep imaging of the field of GRB 10 1225 A in g' and r' 
with OSIRIS at the 10.4m GTC under good conditions and seeing around 1.0 arcsec. The total exposure time is 
8400 s in g' and 6400 s in r' band. We detect an unextended object at the position of the GRB in both bands at 
g’ — 27.21 ± 0.27 (~3 a) and r' = 26.90 ± 0.14 (~la). These values are clearly above the extrapolation of the SN 
light curve to 180 days and we therefore propose this object to be the host galaxy of GRB. 101 225 A. At a redshift 
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Figure 6: Pre-imaging exposure in g' band obtained with the 3.5m CFHT. The field of view is 60 arcsec x 
40 arcsec, North is to the top and East to the left of the image. The circle indicates the position of the optical 
counterpart of GRB 101225 A, where we see a la detection of what could be the host galaxy. 


of z = 0.33, the absolute luminosity is M ab s, g — -13.7 mag, ~ 2 mag fainter than the faintest GRB host detected 
so far (XRF 060218 [12]) and one of the faintest galaxies ever detected at that redshift. The g' - r' color of the 
host is consistent with what we would expect for a late-type star-forming galaxy, in line with what is seen for other 
long GRBs, although we note that the large errors allow other galaxy types as well. The blue color would also 
argue against the source being a late detection of a cool component from the event itself (see [3]). 



Figure 7: Imaging of the host using OSIRIS/GTC 6 months after the GRB. Panel a shows the stacked image 
in /, panel b the same in the g' band. The field of view is 50 arcsec x 65 arcsec. The circle indicates the position 
of the host galaxy. 
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4 Optical spectroscopy 


We obtained a spectrum of the optical counterpart 5 1 h after the event using OSIRIS on the 10.4m GTC telescope 
on La Palma (Spain). Two spectra of 1800 s exposure time each were taken with grism 300B (R=325, wavelength 
range: 3500 - 7000 A) under moderate to high airmass (1.26 and 2.05, respectively). The spectra were reduced 
and combined with standard tasks in IRAF and flux-calibrated with the spectrophotometric standard G191-2B2 
observed in the same night. The continuum is clearly detected, but the spectrum shows no obvious absorption or 
emission lines. 

We performed a search for any possible redshift solution in the range z— 0.1 - 0.6. For this, we stacked the 
normalized spectrum at the position of the strongest emission lines [Oil] AA 3727, 3729 A, [OHI] AA 4959, 5007 
A, H/3 and H a at redshifts between 0.1 and 0.6 in steps of 0.005. We find no indication of any emission feature in 
the stacked spectrum at any redshift in this interval. 

The limits on the detection of Ha [OIH] and [OH] emission from the host galaxy are < 5 x 10~ 18 erg cm -2 s _1 , 
< 2.3xl0 -18 erg cm -2 s -1 and < 3 x 10“ 18 erg cm -2 s -1 (3cr) respectively. We can also put a limit on the de- 
tection of Ha at 2 = 0 of < 2x 10 -18 erg cm -2 s _1 . The flux-calibrated spectrum with the position of typical 
emission lines from the host shifted to a redshift of z = 0.33 is shown in Fig. 8. 



4000 5000 6000 7000 8000 9000 10000 

observed wavelength in A 

Figure 8: Flux-calibrated spectrum of the optical transient from GTC/OSIRIS 2.1 days after the GRB. The 

error spectrum is plotted in blue. The red lines indicate the position of normally strong emission lines from the 
interstellar medium at a redshift of z = 0.33, none of the lines is detected in our spectrum. 

On the night of 201 1 Feb. 04 we observed the optical counterpart with the Low-Resolution Imaging Spectrom- 
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eter [35] on the Keck I telescope during local twilight. Two undithered observations of 600 seconds each were 
acquired using the 1.0 arcsec slit and the D500 dichroic at a position angle of 86.5 degrees. On the red arm we 
used the 600/7500 grating and binned the CCD along the spatial direction (2x1 binning); on the blue arm we used 
the 600/4000 grism and binned the data along both spatial and spectral axes (2x2 binning). Due to twilight there is 
no evidence of a trace in the blue spectrum (and no source is detected in GMOS g'-band imaging from the previous 
night). A faint continuum trace at the expected position of the transient is identified on the red side in the second 
(less twilight-affected) exposure from 7160 to 8000 A with no absorption or emission features visible. 

5 Modeling the UV to NIR spectral energy distribution 


The early evolution of the UV/optical/IR (UVOIR) counterpart is very unusual for a GRB afterglow. Instead of a 
power-law spectrum with a negative spectral slope, it had a very blue counterpart, following what seemed to be a 
power-law with a positive spectral slope [36]. Furthermore, the counterpart stayed bright during the first days and 
then decayed [37] with a strong color change, transforming into a very red counterpart two weeks after the trigger 
[38]. We interpret this early evolution as being produced by the expansion and cooling of a blackbody (BB), as 
shown in Section 5.1. 

The simple BB evolution is not valid any more for the emission beyond ~ 20 days after the trigger. At that 
time we observe a flattening of the light curve, while the very red color is preserved. This late evolution can be 
well -described with the presence of a supernova component, as described in Section 5.2. 

5.1 Early time evolution 


For the modeling of the UVOIR spectral energy distribution (SED), we use the photometry presented in Tab. 2 
together with some of the data points extracted from the literature [36, 37, 38, 39, 40, 41, 42]. All magnitudes are 
corrected for a Galactic extinction of Ay = 0.33 mag and transformed from magnitudes to flux densities. With 
this data set we are able to derive a set of 12 SEDs ranging from 0.07 to 40 days after the trigger. 

The early optical SEDs are well fitted by using an expanding and cooling blackbody of the following form (in 
frequency space): 


f„(Jv) - io 26 


f R\ 2 2-irhv 3 (1 + z) 4 1 


( 1 ) 
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Here, the factor 10 26 is used to convert W/m 2 /Hz to Jy. R is the radius of the emitting black body (which we 
assume to be spherical), D is the luminosity distance to the object, 2 the redshift and T 0 bs is the observed blackbody 
temperature (the rest-frame temperature would be T rest = T obe { 1 4- z)). The other physical constants are: c is the 
speed of light, h Planck’s constant and ks Boltzmann’s constant. For simplicity we assume a blackbody with an 
emissivity of 1 . The blackbody succeeds in reproducing the data up to 10 days, without any intrinsic extinction or 
additional emission component, after which another component becomes dominant (see Section 5.2), 


Table 4: Measured values for the blackbody evolution. Values in brackets are extrapolated from the BB evolution 
due to a limited amount of data points in those SEDs. 


Epoch 

(days) 

Observed temperature 
(K) 

Normalisation constant 

(1 + Z )W«(#) 2 

0.07 

(43 000±8 000 

1.7±1.5) 

0.17 

40 000±6 000 

1.8±1.2 

0.3 

35 000±3 950 

2.2±0.8 

0.6 

25 340±5 440 

4.6±2.5 

1.1 

20 900±1 770 

5.2±1.2 

2.0 

15 000± 1 090 

8.9±1.8 

3.0 

14 260±1 760 

8.00±3.4 

5.0 

(11 300±2 000 

10.0±4.0) 

10.0 

(6 000±2 000 

14.0±4.0) 

18.0 

(5 000±1 000 

18.04= 10.0) 


From the fits to the SED evolution and allowing for a second-order fit, we get the following evolution of the 
normalisation constant: 

log ^10 26 tt(1 + z ) 4 ^ = (0-70 ± 0.04) + (0.46 ± 0.03)Joff(i) - (0.01 ± 0.0b)logl0(t) 2 (2) 

where t is the time in days. The temperature evolution (in K) can be described by: 

log (T oba ) - (4.342 ± 0.017) - (0.395 ± 0.016 )log{t) - (0.11 ± 0.02)log{t) 2 (3) 

Figure 9 shows the temporal evolution of the normalisation constant and the temperature. The normalisation 
can be sufficiently described by a linear evolution in log-log space and therefore the second order term in eq. (2) 
can be neglected. For the temperature, we need an additional second order term to obtain a reasonable fit to the 
data. The temperature and normalization evolution fits very well to our theoretical model of the afterglow (see Sect. 
3 and Fig. 19). Once the redshift is known the normalisation constant can be transformed into physical values in 
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the rest-frame of the object which is shown in Fig. 10. For this we assume z = 0.33 (see Section 5.2), or 1661.1 
Mpc using a ACDM cosmology with Hq = 71, f Im - 0.27 and Qa - 0.73. 



0.1 1.0 10.0 
observed time (days) 


Figure 9: Temporal evolution of the normalization constant and temperature during the UVOIR BB phase. 

The evolution of the normalization constant can be described by a linear evolution in log-log space while the 
temperature evolution requires a second order term. 


5.2 Late evolution and SN template fitting 

The late evolution of the light curve requires a component in addition to the evolving BB described in Section 
5.1. A pure BB fit to the late data gives a bad fit with x^/d.o.f. - 10.45/5 and is therefore ruled out with 
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Figure 10: Temporal evolution of the physical radius and temperature of X-ray and UVOIR BB. This figure 
is identical to 9 but showing the evolution of both BB components together, 

a 94% probability. We assume the late SED to have an additional component from a SN and fit the SED with 
several SN templates. These fits are also used to estimate the redshift of GRB 10 1225 A, which we could not obtain 
spectroscopically (Section 4). 

To determine the redshift of GRB 101225A we use the SED at 40 days after the burst where the contamination 
from the BB is negligible and where we have detections in 7 different bands. Given the steep slope in the blue part 
of the spectrum, we convolve the response of each filter with the spectral shape of the templates. This is particularly 
important for the r'-band observations performed from GTC and Gemini at a very similar epoch, which show a 
significant difference in flux density. The filter of GTC reaches slightly redder wavelengths, and the difference in 
flux densities can be well explained by a very steep slope due to a SN feature as shown in Fig. 12. 

We obtain templates for different core-collapse supemovae from the literature 1 . We exclude SN la from the 

' http://supemova.lbl.gov/~nugent/nugent_templates.htiii! 
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Figure 1 1 : Evolution of the velocity of the UVOIR black body according to the result of our modeling. The 

plot shows the evolution at the best-fit redshift of z = 0.33^ 0^4 (red dashed line and grey dotted lines) as well as 
three other redshifts. At z = 0.33, the velocities derived during the BB phase range between 0.25 - 0.005c, typical 
of accelerated material during a supernova explosion. 


analysis, as we do not expect high-energy emission such as detected for GRB 101225A for those events. The 
template for each SN prototype was interpolated to the time of the SED for a range of redshifts (see references 
in Table 5). In the particular case of SN 1998S the templates were created by combining ground-based and HST 
spectra. 

In order to obtain the best-fit redshift, we also need to consider the time evolution and maximum of the SN 
light curve which is expressed as the stretching factor (a stretching factor of s = 1 corresponds to a time evolution 
identical to the corresponding template SN). To this end, we iteratively fit the SED and the light curve starting with 
s = 1 and a first fit to the SED in 7 bands (see Fig. 12). With the best-fit z , a SN light curve is derived and fitted 
to the real light curve which gives a new stretching factor. This process continues until the solutions for s and z 
converge. We also allow for a scaling of the flux of the SN during the fit. The best fits for s and z for each SN 
template are displayed in Table 5. 

The absolute best fit is obtained with a SN 1 998bw template, a broad-lined Type Ic that is the canonical template 
for GRB-related supemovae. For this case we obtain a redshift of 2 = 0.33 + q;qJ and s = 1.25 ± 0.15. The other 
core-collapse SNe we tested (excluding the Type II SN 1998S, which clearly does not fit our SED) give redshifts 
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Table 5: Fits of the SED with SN templates. A pure BB fit (not listed in this table) is rejected with a probability of 


94%. 


SN Type 

SN template 

best fit z 

stretching factor 

X z /d.o.f. 

rejection prob. 

reference 

Ib/c 

1999ex 

0 3i+u.ut) 

1.1+0.15 

5.80/5 

67% 

[43] 

Ic broad-lined 

1998bw 

0.33lS;SI 

1.25+0.15 

3.60/5 

39% 

[44,45] 

n 

1998S 

O.501SS 

1.0+0.15 

44.5/5 

100% 

[46, 47, 48] 

iin 

200 laj 

— 

— 

> 50.0/5 

100% 

[49] 

EL 

1985P 

0 AT. +0.04 

1.3 ±0.15 

6.1/5 

70% 

[50] 

IIP 

1999em 

0 41+0-05 
u/fi -0.03 

1.0 ±0.15 

6.3/5 

72% 

[51] 


between z = 0.31 and z = 0.50. We therefore use z — 0.33 as reference in this work. Figure 12 shows the fitted 
SED with the template at z - 0.331q ol an d Figure 13 the best fits for all other SN templates as well as a simple 
BB. 



2x1 0 14 3x10 14 5x10” 8x10” 

frequency (Hz) 


Figure 12: Fit of the day 40 SED to a SN 1998bw template. The observations are in black and the best fit, with 
a redshift of z - 0.33, in red. The gray lines represent the template at the different redshifts (steps of 0.01) within 
errors. 


5.3 Luminosity and stretching factor of the SN associated with GRB 101225A in context 
of other GRB-SNe 


We undertake a more general comparison to SNe associated with GRBs by following the formalism of [52]. 
These authors used a SN 1998bw template light curve to fit late bumps in GRB afterglow light curves, modifying 
the template by increasing or decreasing the luminosity at peak (the parameter k, with k 1 implying a peak 
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frequency (Hz) 

Figure 13: Fit of the SED at 40 days to SN templates other than a SN 1998bw-like (as listed in Tab. 5) The fit 
shown is the one with the best-fit stretching factor for each template. We also plot a simple BB with 3650 K, the 
temperature at which we obtain the best fit. 

luminosity identical to that of SN 1998bw), and stretching or compressing the light curve in time while retaining 
the overall shape (the parameter s, again, s = 1 implies the temporal evolution is identical to that of SN 1998bw 
in the same band). This procedure also included the creation of synthetic templates by interpolating between the 
SN 1998bw light curve in different filters, and taking into account the cosmological it -correction. Nearly all 
GRB-SNe were well-fit by the SN 1998bw light curve template. For GRB 101225A we have s = 1.25 ± 0.15 and 
k = 0.08 ± 0.03 according to the designation of [52] using the SN 1998bw light curve at z = 0.33. 

Ferrero et al. [53] analyzed SN 2006aj associated with XRF 060218, and placed it into the k — s context. They 
employed the line-of-sight extinction values derived by [54] to derive intrinsic k values. To place the SN associated 
with GRB 10 1225 A into the k — s context, we fit the light curve analogous to [52], and use the sample of [53] as 
well as additional events as a comparison. The complete data are presented in Table 6. 

GRB 990712 has been analysed again with additional data. We find no evidence for host extinction. For GRB 
02121 1, a re-analysis of the afterglow SED finds no evidence for host extinction, the value from [53] thus remains 
unchanged but now counts as extinction-corrected. For GRB 040924, we use the k and s values from [9] and 
correct k with the extinction found by [54]. For GRB 050525 A, we use the uncorrected k value from [53], and 
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correct it with the extinction found by [11]. For XRF 050824, we use the uncorrected k value from [55], and 
correct it with the extinction found by [11]. GRB 060729 is analysed in [56], GRB 0803 19B is analysed in [57]. 
GRB 090618 has been analysed for this work, using the data set of [18]. We were not able to derive a good SED 
for this afterglow, therefore the k value has not been corrected. 

As can be seen in Fig. 14, the SN associated with GRB 101225A is significantly fainter than any other known 
GRB SN (with the SN associated with GRB 040924 being the most similar, but this event is only marginally 
detected). At the same time its temporal evolution is similar but slightly slower than most known GRB-SNe, 
though not by a large amount. It is also fainter than two well-studied Type Ic SNe, SN 19941 and SN 2002ap (see 
[53] for discussion), the latter being broad-lined, but not associated with a GRB. 
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Figure 14: Luminosity factor k and stretching factor s of SNe associated with GRBs. Filled symbols have 
been corrected for host-galaxy line-of-sight extinction, non-filled symbols have not. We label several well-studied 
nearby GRB-SNe, as well as two “canonical” Type Ic SNe, SN 19941 (Ic) and SN 2002ap (broad-lined Ic unasso- 
ciated with a GRB). The GRB 101225A SN is fainter than all these events. 


6 Discussion on the redshift of GRB 101225 A 


Determining the distance scale at which GRB 101225A occurred is crucial to understand the energetics and get a 
clear picture of the physics involved in this event. In this section, we present several independent arguments to 
strengthen our redshift estimation. 

The first strong limit on the redshift comes from the UVOT detection in uvw 2 which implies a redshift lower 
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Table 6: Luminosity Factor k and Stretching Factor s for GRB SNe 


GRB 

k 

s 

Comment 

Reference 

970228 

0.41 ± 0.29 

1.45 ±0.93 

uncorrected 

[53] 

990712 

0.60 ±0.10 

0.76 ± 0.07 


This Work 

000911 

0.85lgS 

1.40 ±0.32 


[53] 

010921 

1 fic:+2.82 

l.OO_ 0 79 

0.69 ± 0.25 


[53] 

011121 


0.30 ±0.02 


[53] 

020405 

O.90±°;i5 

0.97 ± 0.07 


[53] 

020903 

0.62 ± 0.09 

0.92 ± 0.08 

uncorrected 

[53] 

021211 

0.43 ± 0.21 

0.99 ± 0.25 


[53] , This Work 

030329 

i p;n +0 - 19 
i..OU_o,16 

0.85 ±0.10 


[53] 

031203 

i oQ+o.is 

j— ^o_0.16 

1.09 ±0.07 


[53] 

040924 

0.24t°;g 

1.37 ±0.97 


[10], [54] 

041006 

l-M+Sg 

1.38 ±0.06 


[53] 

050525A 

0.68iS;g 

0.77 ± 0.04 


[53], [11] 

050824 

i 20 + °' 69 
i.zu_ 0 56 

0.52 ±0.14 


[55], [11] 

060218 

0.74 ±0.01 

0.69 ±0.01 


[53] 

060729 

1.04 ±0.03 

0.86 ±0.02 


[56] 

0803 19B 

Z.OU— 0 90 

0.89 ± 0.10 


[57] 

090618 

0.58 ± 0.05 

0.99 ± 0.05 

uncorrected 

This Work, [19] 

101225A 

0.08 ± 0.03 

1.25 ±0.15 


This Work 


than z - 1.4 [58]. As an independent test, the analysis of the absorption in the X-ray spectra (see Sect. 2 and 
in particular Fig. 4) imposes an upper limit on the redshift of 0.5 within 99% confidence and 0.35 within 90% 
confidence. 

The most restrictive redshift estimate comes from the SN fitting at 40 .days after the burst (see Sect. 5.2), 
where we find a best fit for a broad-lined SN Ic, such as SN 1998bw, the prototype of GRB-related SNe, at 
2 = 0.33±g;gJ. For the rest of the SN types that give a reasonable fit, we obtain similar values, always resulting in 
a redshift between 0.21 and 0.50. Independently of the SN fit, we can compare the brightness of the bump in the 
light curve with the dimmest known SN (SN2008ha, which peaked at M R = -14.5 ± 0.3, [59]) and take that as a 
lower limit on the redshift. In this case, we would expect GRB 101 225 A to be located at a redshift larger than 0.1. 

From the SED fit of the first days, we know that the evolution is well-described by a simple BB. Depending on 
the distance at which the object is found, we can derive different radii and expansion velocities. For an explosion 
of this type, we expect expansion velocities larger than ~ 10 3 km s _1 , which would be at the edge of a stellar 
wind regime and, if similar to a SN explosion, of the order of 10 4 km s -1 . We cannot, in principle, rule out 
higher velocities of the ejecta. Under the assumption that the ejecta should not be traveling at velocities larger 
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than 100,000 km s -1 (0.3c) and lower than 1,000 km s" 1 (0.03c), we can estimate a range of reasonable velocities 
between z ~ 0.20 and z ~ 0.60 (see Fig. 1 1). We note that there might be additional effects such as variations in 
the transparency of the BB that could introduce variations in this simplified analysis. In any case, the evolution of 
the BB would be hard to get by any source within the Local Group. At a redshift of z = 0.33, the velocity of the 
blackbody would have evolved from ~ 70, 000 km s _1 at the time of our first SED to ~ 2, 000 km s -1 at the time 
of our last blackbody-dominated epoch, nicely matching the requirements. The fact that the blackbody detected for 
XRF 06021 8/SN 2006aj was very similar to the one found for GRB 101225Awhen placed at a redshift of z = 0.33 
(see Section 7) adds additional evidence for the validity of the redshift estimate. 

Concluding, the proposed redshift for GRB 101225A of 2 = 0. 33^^04 * s supported by several independent 
arguments, and can be considered as a firm reference when studying the physical processes involved in the event. 
We do not find any evidence that would indicate a redshift smaller than 2 = 0.2 or higher than 2 = 0.5. 

7 Comparison between GRB 101225A and other GRBs with SNe and BB 
components 

In Table 7, we compare several properties of GRBs associated with SNe and without a “classical” afterglow com- 
ponent. All of them are subluminous compared to the average long-duration GRB with Ei S0 around 10 51 — 10 54 
erg. GRB 101225A lies on the lower end of the energy output from normal long-duration GRBs. Among those 
nearby GRB-SNe without a classical afterglow, there is a class of very long duration GRBs with very low E pea ^ 
values, all of them showing a thermal component in X-rays. XRF 060218 and XRO 080109 also had a thermal 
component at optical wavelengths during the first few days [13, 22]. For XRF 100316D, no optical counterpart 
was detected before the onset of the actual SN due to high intrinsic extinction in the host galaxy. 

XRF 060218 [13] shows a similar early behavior to GRB 101 225 A. We compare the early UVOT light curve of 
XRF 060218, which we obtained from the UVOT catalogue [60], to the light curve from GRB 101225A by shifting 
that of XRF 060218 to 2 = 0.3, including a AT -correction. To obtain the K -correction, we use XSPEC assuming 
a blackbody spectrum with kT ~ 3.7 eV. The temperature was determined from the best-fit model of a BB to a 
SED of XRF 060218 taken at 120ks after the trigger [13]. Using this BB spectrum, we determined the expected 
flux density in the observed frame for each filter and at 2 = 0.33. The ratio of these two flux densities was taken to 
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Table 7: GRBs with SNe but without afterglows 


GRB 

z 

TC? 

T 90 

(s) 

E peak 

(keV) 

Ejso 

(erg) 

HR(50-100) 

Radio? 

(mag) 

SN My 
(mag) 

Host Mb 

980425 

0.0085 

No 

23.3 

55±21 

8.1 xlO 47 

— 

Yes 

-19.42 

-17.6 

031203 

0.105 

No 

30 

158±51 

3xl0 49 

— 

Yes 

-20.39 

- 21.0 

060218 

0.0331 

Yes 

-2100 

4.9 

6 . 2 xl 0 49 

0.835 

Yes 

-18.76 

-15.9 

080109* 

0.0065 

Yes 

—400 

low 

2 xl 0 46 


Yes 

-16.7 

-20.7 

100316D 

0.059 

Yes 

>1300 

— 

3.1 xlO 49 

0.891 

No 

-18.62 

-18.8 

101 225 A 

0.33 

Yes 

>2000 

38±20 

> 1.4 x 10 51 

1.06 

No 

-16.9 

-13.7 


No 7 -rays observed, numbers derived from X-rays. 

TC? refers to the early thermal component, mostly attributed to a supernova breakout. 
HR is the hardness ratio, defined as the ratio of channels (50 — 100 kev)/(25 — 50 keV) 
SN My is the SN peak absolute magnitude in V 
Host Mjg is the host absolute magnitude in B 


be the i^-correction for the specific filter, which is — 2.20 for all filters. Figure 15 compares the flux-density light 
curves of the 3 UV filters of both GRBs. For both GRBs the light curves were corrected for Galactic extinction. 



Figure 15: UV light curves in flux density for GRB 101225 A and XRF 060218. The light curves for XRF 
060218 have been shifted to z = 0.33 for direct comparison. The colored solid shapes connected by solid lines are 
for GRB 101225A, open symbols and dotted lines for XRF 06021 8. Circles are uvwl, squares are uvm2, triangles 
are uvw2. 


We furthermore performed a fit of the early UVOIR SEDs of XRF 060218 and SN 2008D, in a similar way to 
what we did for GRB 101225A. XRF 060218 also seems to follow a blackbody evolution at early times. However, 
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the SN starts to dominate already around 3 days after the burst (see Fig. 16), limiting how long the study of the 
evolution is possible. As can be seen from Fig. 17, the evolution of the blackbody is not very different from what 
we see in GRB 101 225 A, although the radius expansion is slightly steeper. SN 2008D does not have a measurable 
thermal component in the early X-ray data. In the optical, the cooling BB dominates until about 4 days after the 
event when the onset of the SN was observed (see Fig. 16). The temperature of the UVOIR BB of SN 2008D is 
lower (~ 30 000 K) than those of XRF 060218 and GRB 101225A. The radius evolution is considerably steeper 
than for those two events (Fig. 17). The UVOIR BB emission of SN 2008D is therefore likely due to the cooling 
of the initial shock breakout. 

8 Temperature evolution and geometry of the UVOIR BB 

The interaction of an ultra-relativistic, well-collimated jet (having a half-opening angle 9j t o ~ 10°) with the 
previously ejected broad torus can explain the power-law and thermal component in the X-rays of GRB 10 1225 A. 
In Fig. 18 we show four different stages in the evolution of the jet and its interaction with the CE shell. 

Taking the results from CE-simulations, we assume that the density of the common envelope ejecta in a narrow 
funnel around the symmetry axis of the system is much lower than elsewhere in the CE-ejecta. This funnel has 
an opening half- angle 9 / ~ 2°, which permits the passage of the 7-radiation generated by internal shocks in the 
ultra-relativistic jet. Since 6 f < 9j , 0 , most of the jet beam hits the inner boundary of the CE-ejecta at a distance of 
RcE,in — 2.5 x 10 12 cm (panel a of Fig. 18), while only a small fraction of the beam propagates through the funnel 
until it reaches the outer radial boundary of the CE-ejecta (Re e, out — 2.1 x 10 14 cm, panel b). As the central 
jet spine progresses through the ejecta funnel, it interacts with the lateral walls, giving rise to mass entrainment in 
the jet. The additional baryon load and the shear with the funnel walls decelerates the ultra-relativistic jet spine 
extremely quickly and, therefore, no regular afterglow signature is produced. Furthermore, most of the jet beam 
hits regions of the CE-ejecta which are much denser than the central funnel and, thus, its propagation across the 
ejecta is much slower than that of the central spine. As the jet decelerates to moderately relativistic speeds it 
expands laterally. 

The X-ray emission is attributed to the shocks (forward and reverse) produced as the ultra-relativistic jet im- 


34 



II 




Figure 16: SED fit of the early UVOIR data of XRF 060218/SN 2006aj (top) and SN 2008D (bottom). The last 
epoch in both plots already begin to show a strong contribution of the SN. These figures can be directly compared 
to Fig. 1 of the main paper. 

pinges against the inner radial boundary of the CE-ejecta. In our model the cross-sectional radius of the CE-ejecta 
funnel (Rf — RcE,in sin Of ~ 10 11 cm) sets the almost-constant size of the observed X-ray source (blue dots in 
Fig. 10). 

An approximate evolution of the jet spine as it crosses the ejecta funnel can be obtained by applying the model 
of [63], which is based on solving a system of four ordinary differential equations (as a function of the observer’s 
time t) for the Lorentz factor, the rest-mass, the radius and the jet half-opening angle, respectively: 
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Figure 17: Temporal evolution of the radius and temperature for the pre-SN epochs of XR1 06021 8/S N 
2006aj (top) and SN 2008D (bottom). The dotted line shows a polynomial tit u< rhe lemperature and radius 
evolution. This figure can be directly compared to Fig. 10. 
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Figure 18: Cartoon of the different stages in our model of GRB 101225 A. The different stages are explained in 
the text. 


dT dm i-r 2 

dt dt Mjet + em + 2(1 — e)Tm 

dm d/2 _ 2/ , 

— - = — 2irR 2 (1 — cos Q)n(R) 

dt dt 

™ = 0cr( r + Vr 2 -i) 

de c g {r+y/r*=i) 

dt R 


(4) 

(5) 

( 6 ) 
(7) 


where n(R) is the particle number density, P = VT= r - 2 is the velocity, e is the shock-generated thermal 
energy (in the comoving frame) that is radiated, and c s — yj (4T + 1)(T 2 - l)(3r(4I' 2 - l)) -1 is the sound 
speed. The temperature of the jet spine is derived by assuming that the pressure is dominated by the photon 
radiation, so that P(t ) aT 4 (f)/3, with the radiation constant a — 7.56 x 10 -15 ergcm~ s K -4 , and noting that 
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P — (r 2 - l)/(3r)nm p c 2 , (m p being the proton mass) according to [63]. 

For the jet spine propagating through the CE-ejecta funnel we may take the following initial conditions, corre- 
sponding to an ultra-relativistic and cold jet: r* n ~ 100, Mi n ~ 8 x 1O -9 M 0 , Ri n = Rce, in = 2.5 x 10 12 cm, 
and 0i n = 9f = 2°. A moderate radiative efficiency e ~ 0.4 is assumed. The initial jet mass is estimated from 
the observed lower limit of E lt i SO and the initial jet half-opening angle as Mj et = Ej et /(({T — 1)(1 — e) 4- e)c r 2 ), 
where Ej et — E 1<iBO (l — cos Of) ~ 8.5 x 10 47 erg. We assume that the CE-ejecta have a declining particle number 
density of the form n(R) — ni n (RcE,in/R) 2 > with = 4.9 x 10 12 cm -3 . After t ~ 0.05 days the jet spine 
reaches Rce, out at a velocity of v jet, out ~ 0.25c, an opening half-angle of 0 jet, out — 70°, and having plowed 
~ 1.5 x lO -5 Af 0 of the CE-ejecta. 

The temperature and radius evolution of the UVOIR BB (panel c in Fig. 18) can be modeled as the result of 
the further deceleration of the sub-relativistic jet emerging trough Rce, out- Again, we use equations 4-7, but now 
considering a uniform medium with a particle density of n ext ~ 1.6 x 10 9 cm -3 and a higher radiative efficiency 
of e = 1. As initial conditions we take the terminal values of the velocity, mass, radius and opening half-angle of 
the previous evolution through the CE-ejecta funnel. The result can be seen in Fig. 19 (after t ~ 0.05 days). Small 
variations of the initial variables can also roughly fit both the radius and the temperature evolution. We note that, 
even though the theoretical model does not predict an exact power-law for the radius evolution, the deviations from 
the observed R(t) oc t 0,22 are rather small. We also point out that the temperature evolution is not compatible with 
a single power-law as the observational data suggests. 

Finally, we speculate on the fate of the jet wings, i.e., the fraction of the jet that does not cross the CE-ejecta 
funnel. The energy contained in these wings, ~ 10 49 erg, is transferred to the CE-ejecta, and can break out 
through Rce, out on scales of ~few days, almost isotropically, but at relatively low temperatures compared with 
the emerging jet spine. Hence, its observational signature will be probably hidden by the emerging SN light curve 
(panel d) in Fig. 18. 
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Figure 19: Radius (panel a) and temperature (panel b) evolution of a sub-relativistic jet. The initial evolution 
ip to ~ 0.05 days (grey shaded area), is computed with a modified version of the [63] method to incorporate a shell 
profile in which the number density decays with radius as n = rii n (RcE,in / ft) 2 , with n in = 4.9 x 10 12 cm -3 . For 
this part of the evolution, the initial data are: R in — RcE,in = 2.5 x 10 12 cm, E jet — 8.5 x 10 47 erg, ~ 100, 
hn = 0/ = 2°, and e ~ 0.4. 
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